UWThPh-2006-28 



A three-parameter model for the neutrino mass matrix 

W. Grimus^ 1 )* and L. Lavoura^ 

Fakultat fiir Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Wien, Austria 

Universidade Tecnica de Lisboa and Centro de Ffsica Teorica de Partfculas 
Instituto Superior Tecnico, 1049-001 Lisboa, Portugal 

24 May 2007 

Abstract 

Using the type-II seesaw mechanism with three Higgs doublets <p a (a = e, fj,, r) 
and four Higgs triplets, we build a model for lepton mixing based on a 384-element 
horizontal symmetry group, generated by the permutation group £3 and by six Z2 
transformations. The charged-lepton mass matrix is diagonal; the symmetries of 
the model would require all the three masses m a to be equal, but different vacuum 
expectation values of the 4> a allow the m a to split. The number of parameters in 
the Majorana neutrino mass matrix M v depends on two options: full breaking of 
the permutation group S3, or leaving a [i-t interchange symmetry intact; and hard 
or spontaneous violation of CP. We discuss in detail the case with the minimal 
number of three parameters, wherein M. v is real, symmetric under \i-r interchange, 
and has equal diagonal elements. In that case, CP is conserved in lepton mixing, 
atmospheric neutrino mixing is maximal, and #13 = 0; moreover, the type of neutrino 
mass spectrum and the absolute neutrino mass scale are sensitive functions of the 
solar mixing angle. 
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1 Introduction 



There are two puzzles associated with neutrinos: why are their masses so much smaller 
than those of the charged fermions, and why does the lepton mixing matrix feature large 
mixing angles — for reviews see P] — in contrast to the quark mixing matrix. It is possible 
that both puzzles are solved through the same mechanism. In this paper we envisage 
the type-II seesaw mechanism as a possible solutionis We use horizontal symmetries to 
enforce certain features of lepton mixing, in particular maximal atmospheric neutrino 
mixing and #13 = 0. In order to achieve this, we enlarge the scalar sector of the Standard 
Model by adding to it four Higgs triplets and by using altogether three Higgs doublets. 
Our model has a permutation group S3 together with six cyclic symmetries Z2, which 
commute with each other but not with 53; the result is a large discrete symmetry group 
with 384 elements. This setting allows to obtain four different neutrino mass matrices, 
depending on the assumed breaking of the horizontal symmetries and of the symmetry 
CP. Amazingly, by breaking the horizontal symmetries softly by terms of dimension 
two, while leaving a residual /x-r interchange symmetry to be broken at low energy in 
the charged-lepton sector, we arrive at a viable neutrino mass matrix with only three real 
parameters. 

In section [2] we make a general discussion of the type-II seesaw mechanism with an 
arbitrary number of Higgs doublets and triplets. Our model, with its multiplets, symme- 
tries, and Lagrangian is explained in section [3j In section H] we investigate in detail the 
most predictive case of a three-parameter neutrino mass matrix. A generalization thererof 
is considred in section [51 The conclusions are presented in section [61 



2 The type-II seesaw mechanism 

We first review the type-II seesaw mechanism [21 [3] for small neutrino masses. We assume 
the existence — in the electroweak theory — of several Higgs doublets <p a with hypercharge 
1/2, together with several Higgs triplets A, with hypercharge 1. Let the neutral compo- 
nents of the 4> a have vacuum expectation values (VEVs) v a and the neutral components 
of the Aj have VEVs Si. Just because of the hypercharge symmetry, the vacuum potential 
Vo must be of the forrru 

V ° = W) a/3 V ° V P + (^a) ij S i S 3 + (UapS*V a Vp + CC.) 

+Kp~f&v* a vpv*vs + \i jk iS*SjSlSi + KpijvlvpSlSj. (1) 

The matrices and fi\ are Hermitian and, likewise, the A coefficients must obey various 
conditions in order that Vq should be real. The VEVs of the triplets are determined by 

= = + tiapVaVp + ^ijklSjS* k Si + \ a f3ijV* a V pS j . (2) 



1 In our model we do not allow for a type-I seesaw mechanism; we assume right-handed neutrino 
singlets not to exist. 

2 We use the summation convention. 
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Contrary to we assume the matrix [i\ to be positive definite so that, in the absence 
of the t ia/ 3 terms, the only solution to equations ([2]) would be for all the 5,i to vanish. The 
VEVs v a are of order of the electroweak scale v ~ 174 GeV, or smaller. Assuming the t ia p 
to be of order M and the eigenvalues of n\ to be of order M 2 , where M is a mass scale 
much larger than v [3], the approximate solution to equations ([2]) is given by [5] 

Si ~ - tjapVaVp. (3) 

From equation ([3]), the 5i are of order v 2 /M <C v. If, furthermore, all the A coefficients 
are of order unity or smaller, then the approximate solution ([3]) will be corrected on its 
right-hand side only by terms suppressed by a factor v 2 /M 2 <C 1. 

Under an 577(2) gauge transformation, the left-handed lepton doublets D^ a trans- 
form as Dlo — > WD La while the Higgs triplets transform as A, — > H^AjH^, where W 
is an SU(2) matrix. Therefore, the Higgs triplets have Yukawa couplings of the form 
D^C^eAiDip, where C is the charge-conjugation matrix in Dirac space and e is the 
2x2 antisymmetric matrix in gauge-S77(2) space. The VEVs Si being very small, the 
above Yukawa couplings generate very small neutrino mass terms <5j ^L a C~ x vi,p, of order 
v 2 /M times a typical Yukawa-coupling constant. The neutrino masses being of order 
0.1 eV, M could easily be of order 10 14 GeV [3], thus fully justifying the approximate 
solution (j3J). 

3 The model 

Our model follows closely, in the symmetries that it utilizes, a previous model of ours [B] . 
We have three left-handed lepton doublets D L „, three right-handed charged-lepton singlets 
an, and three Higgs doublets 4> a (a = e, //, r) There are four Higgs triplets, A a and A 4 @ 
The symmetries of the model consist of a permutation group S3 acting simultaneously on 
all indices a, three symmetries 

: 4>a -> ~4>a, «fl -> -®R, (4) 

and another three symmetries 

zi 2) : D La -> -7J L a, «i? -> -or, and A^ -> -A^ iff /3 7^ a. (5) 

Notice that A 4 is invariant under all these symmetries. In appendix A we make a study 
of the full symmetry group of our model. 

The Yukawa Lagrangian invariant under all these symmetries is 

I 

^Yukawa = ~VqD La (j) a a R + - yiDl a C' 1 e A 4 D La 

+y 2 {D T Le C-hA^D LT + D T Lp C- x eA T D Le + D T Lr C- x eA e D Lv ) + H.c. (6) 

3 Constraints on multi-Higgs doublet models from electroweak precision tests are not very stringent: 
Higgs bosons with large ZZ couplings must have an average mass in the range allowed for the mass of 
the Standard Model Higgs boson [4]. 

4 The scalar content of our model resembles that of the A4 model of [7]. However, in that model, three 
gauge triplets are used instead of our A4. 
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Thus, the charged-lepton mass matrix is automatically diagonal, the charged lepton a 
having mass m a = \yoV a \. On the other hand, the neutrino mass matrix is 



yih V-iSr 2/2^ 

M„ I y 2 S T yi5 4 y 2 S e 
y 2 5 il y 2 S e yi$4 



(7) 



all its diagonal matrix elements being equal. 

Due to the symmetries of our model, the coupling constants ti a p of the previous section 
assume the very simple form 



U 



(8) 



Hence, from equation 



Si = -tv a v a (/4) i4 • (9) 
Ordering the triplet fields as (A e , A M , A r , A4), the symmetries of our model would enforce 

(& = dia S (a*1) lAi tfi, t4) > ( 10 ) 

which is not satisfactory since it would lead, through equation (jHJ), to 5 e = 5^ = S T = 0. 
We must have 7^ for i = e,fj,,r. In order to solve this problem, we assume the 

symmetries of the model to be broken softly, only by terms of dimension two. Without 
any residual symmetry, this means that both matrices fx\ and become fully general, 
while all other couplings remain unchanged. 

However, in order to simplify our model and render it more predictive, we may assume 
that the interchange symmetry \l «-> r [8j [9j QUI EH] , which is a subgroup of our permutation 
group S3, is kept unbroken in n\ and Then, 



(A)' 
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(a, d, e, and g are real), so that S e = —tcv a v a , 5^ = 5 T = —tfv a v a , and the neutrino 
mass matrix is /x-r symmetric. This immediately leads to the predictions 623 = 7r/4 
and 6*i3 = 0. The \i-t interchange symmetry is supposed to be spontaneously broken 
through the VEVs of the Higgs doublets: t> M 7^ v T . The Higgs potential is rich enough to 
allow for this outcome — in appendix B we demonstrate this by working out a simplified 
case. Of course, the spontaneous breaking at the electroweak scale of the /z-r interchange 
symmetry will seep, through radiative corrections, into the rest of the theory, so that at 
loop level the matrix (fJ,\) 1 will not any more be of the form in equation (fill , and then 
^ ^ i T . But, both because this is a loop effect, and because it is a correction of order of 
the ratio of the electroweak scale to the much larger mass terms in /i^, we may expect 
Sfj, — S T to remain negligible. 

In a further simplification of our model, we may also assume CP violation to be 
spontaneous: the matrix (ft^) 1 is then real, but the VEVs v a display non-trivial relative 
phases. Then, S e , 5 M , S T , and <5 4 will all have the same phase — the phase of v a v a . That 
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phase may be rephased away from A4 U , so that the neutrino mass matrix becomes real. 
Thus, spontaneous CP breaking in our model yields the remarkable outcome that, even 
though there is CP violation, it remains absent from the mass matrices and from lepton 
mixing!^ 

One thus obtains the following four possibilities: 

1. The general case, in which CP violation is hard and /i-r symmetry is allowed to be 
broken in fi 2 A . Then, 



with real m, p, q, and r. This case should not be very predictive, since it has seven 
parameters to predict nine observables — three neutrino masses, three lepton mixing 
angles, one CKM-type phase, and two Majorana phases. 

2. The case in which /j-t symmetry is allowed to be broken in fi A but CP violation is 
spontaneous. Then, ip, x, and p in ( fT2l) vanish. There is no CP violation in lepton 
mixing. The four parameters m, p, q, and r allow one to predict six observables — 
three neutrino masses and three lepton mixing angles. 

3. The case in which /i-r interchange symmetry is preserved in fi 2 A , while CP violation 
is allowed to be hard. Then, 



with complex parameters x, y, and w. There are in this case five parameters — three 
moduli and two phases. 

4. The most predictive case, in which CP violation is spontaneous and /i-r interchange 
symmetry is preserved down to the electroweak scale. The neutrino mass matrix is 
the one in equation ([13]) but with real x, y, and w. The neutrino mass matrix has 
only three parameters. 

4 The three-parameter neutrino mass matrix 

In this section we concentrate on case 4 of the previous section, i.e. on the neutrino mass 
matrix of equation ( fTBl with real x, y, and w. The algebra of the diagonalization of a 
general /i-r symmetric neutrino mass matrix has been worked out in [9], and we only 
need to adapt it to the simpler case 4. In the following, the solar mixing angle — which 
is defined to be in the first quadrant — is denoted 9, the neutrino masses are m 12 ,3, the 

5 This is not an original situation; in the classical Branco model of CP violation [12], spontaneous CP 
breaking also does not find a way into the quark mixing matrix. 




(12) 




(13) 
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solar mass-squared difference is Am^ = m| — m\ > 0, and the atmospheric mass-squared 
difference is 

eU-^±Hi), (14) 



^ m atm 



2 m? + m 2 , 




where e = +1 indicates a normal neutrino mass ordering and e = — 1 an inverted ordering. 
Equations (3.9)— (3. 11) and (3.15) of [9] yield, respectively, 

m 3 = \x — w\ , (15) 



x 2 + Ay 2 + (x + w) 2 =F Am 



™i, 2 = ~ 2 -, (16) 

tan 20 = 2v ^\ y \ S ign(2x + w), (17) 
Am 2 :) cos26' = w (2x + w) . (18) 



Experimentally we know that 9 is in the first octant. Hence tan 29 > and equations (p3 
and (fl8|) both give 

sign (2x + w) = signw. (19) 
Again from equations ( fT7|) and ( TTSl . 

Iw|tan2# 

M - Ht^- (20) 

Am'i cos 29 — w 2 . . 

x = 2 . 21 

2w v ; 

From equations (Tl4l) . (|T5|) . and (|T6l) . we find 

eAm 2 tm = — - 3xw - 2y 2 . (22) 
Inserting equations fl20l) and (1211 into equation fl22l . we obtain the value of w: 



W " 8 - tan 2 29 " ^ 



4eAm 2 tm + 6Am| cos 26> 
8 - tan 2 29 

Since Am 2 tm ^> Am^, the numerator of equation (|23|) has the sign of e; hence its de- 
nominator must also have the sign of e. That denominator vanishes when sin 2 # = 1/3, 
i.e. when 9 is just the Harrison-Perkins-Scott (HPS) solar mixing angle [10]. We thus 
conclude that, in our model, 

• if the neutrino mass spectrum is normal, then the solar mixing angle is smaller than 
its HPS value; 

• if the neutrino mass spectrum is inverted, then sin 2 9 > 1/3. 

This remarkable result relates the type of neutrino mass spectrum to the value of the 
solar mixing angle. 
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From equations (jT5jl, (fTBTh (T2U]h and (j2Tj). 



1 / |w| AmQCOs2#\ 3|u>| AmQCOs2# 



m lj2 = - f-j and m 3 = . (24) 

2 \cos2ft \w\ J 2 2\w\ 

This gives the overall scale of the neutrino masses. Since \w\ diverges when tan 2 26 — > 8, 
we see that in our model the neutrino mass spectrum becomes quasi- degenerate when the 
solar mixing angle approaches its Harrison-Perkins-Scott value. 

One easily sees the reason why our model displays a singularity when sin 2 9 = 1/3. 
The most general neutrino mass matrix leading to HPS lepton mixing is 

/ x y y \ 
M u = y x + u y — u . (25) 
\ y y — u x + u J 

Our M. v in equation (TIB"]) has equal diagonal matrix elements. Hence, if it were to accept 
sin 2 9 = 1/3, it would have to correspond to u = in equation (125]) . But the M. v of 
equation (1251) with u = leads to two equal neutrino masses, hence it is unrealistic. 
There is therefore a contradiction with experiment in the assumption that our Ai u of 
equation ([TBI might be compatible with HPS mixing. 

Experimentally sin 2 9 is close to 1/3, therefore there is the danger that our neutrino 
masses are too large and saturate the cosmological bound [13]. As a numerical exercise, 
we take the 1 a bound on solar mixing from 



0.27 < sin 2 9 < 0.32 3.73 < tan 2 29 < 6.72. (26) 

The mean value of 9 is given by sin 2 9 = 0.30 and tan 2 29 = 5.25. Note that the upper 
2 a limit sin 2 9 = 0.36 gives tan 2 29 = 11.76, which is already significantly larger than 8. 
Thus, there is experimentally ample room for the neutrino masses to be sufficiently small. 
This happens because tan 2 29 is a rapidly varying function of 9. 

In figured] we have plotted m 1; m.3, and mi + m<i + 777,3 against sin 2 9 in our model. 
We have used the best-fit values Am 2 tm = 2.5 x 10~ 3 eV 2 and Am| = 7.9 x 10 _5 eV 2 
from [14J; for 9 we have used the 3 a bounds 0.24 < sin 2 9 < 0.40, from the same sourceJ§ 

Another important observable is mpp, the effective mass relevant for neutrinoless 2/3 
decay. This is equal to the modulus of the (e, e) matrix element of M. V1 i.e., in our case, 
to \x\. Thus, 

|ti?| Am 2 ?1 cos26 l , , 

m * = V — sjsr— < 27 > 

Since 




Am 



2 



¥ _ a5 L_» V A m |cos2 9 , (28) 

we see that in our model we have the relation 

mpp w m 3 /3. (29) 



6 A different model in which the neutrino mass spectrum is normal or inverted depending on whether 
sin 2 9 is smaller or larger than 1/3, and the neutrinos become degenerate in the limit sin 2 9 — » 1/3, has 
been suggested in [15] . 
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Figure 1: Plot of mi, 777.3, and Y2i rrii as a function of sin 2 9 in our three-parameter model. 
The values of Am 2 tm and Awq have been fixed at 2.5 x 10~ 3 and 7.9 x 10~ 5 , respectively, 
in eV 2 . The dashed vertical line indicates the singularity at sin 2 # = 1/3. 



This same relation has recently been obtained in a different model [16J. 

One may ask oneself whether the neutrino mass matrix displays any characteristic 
texture in the limit tan 2 29 — > 8. A glance at equations f[2"Uj) . (|21[) . and (12"3"1) allows one 
to conclude that, in that limit, all the matrix elements of M. v diverge. Moreover, from 
equations (I2T!) and (I20l) . respectively, we obtain 



W 2' [if I 



1 (30) 
for tan 2 29 — > 8, from where the texture of M. v in that limit can be read off. 



5 Extension to the complex case 

In this section we investigate what happens when one allows the neutrino mass matrix 
of equation ( flBl) to have complex matrix elements. Does the intriguing feature of quasi- 
degenerate neutrinos in the limit of HPS mixing, found in the previous section for the 
case of real matrix elements, still hold true? We shall see that it does not; indeed, the 
general neutrino mass matrix ffTBl does not seem to have much predictive power beyond 
Ue3 = and maximal atmospheric neutrino mixing. 
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The symmetric matrix 

M v = I y z w I (31) 




is diagonalized in the following way: 

U T M V U = diag {m 1 ,m 2 , m 3 ) , (32) 

the matrix U being unitary while the rrij (j = 1, 2, 3) are real and non- negative. Due to 
the special form of Ai u , wherein (A4 U ) 12 = (M. u ) 13 and (A4 U ) 22 = (M u ) 33 , U is of the 
form 

/ c s \ 

C/ = diag (e^, e il3 ) I -rs rc r diag (e iS \ e iEa , e iSa ) , (33) 

V — rs rc — r / 

where c = cos6*, s = sm9, and r = 2 -1 / 2 . From equations (I31~]) - (l33l . we find 

x = e" 2 ^ (c 2 m ie - 2iEl + s 2 m 2 e- 2iS2 ) , (34) 
z = — ^— (s 2 m ie - 2jSl + c 2 m 2 e- 2jS2 + m 3 e _2tS3 ) . (35) 

We define rhj = m.,e _2jSj for j = 1, 2, 3. We also define x = 2 (i? — y). Then, the condition 
x = 2, which makes the .Mt, of equation (I3TI) identical with the one of equation ({TBI) , is 
equivalent to 

mi (2c 2 e ix - s 2 ) + m 2 (2s 2 e ix - c 2 ) - m 3 = 0. (36) 

Thus, the condition x = z is equivalent to the existence of four phases x an d ^1,2,3 such 
that the condition (|36|) is satisfied. That condition states that it is possible to draw a 
triangle in the complex plane, the sides of that triangle having lengths s/Amx, \^Bm 2 , 
and 777.3, where 

A = 4c 4 + s 4 — 4c 2 s 2 cos Xi , x 

(37) 

B = 4s 4 + c 4 - 4c 2 s 2 cos x- 

Therefore, one may eliminate the phases £1,2,3 from condition (1361) by writing the sole 
"triangle inequality" [17] 



A 2 m\ + B 2 m\ + m 4 - 2 (ABm\m\ + Am\m\ + Bm 2 2 mfj < 0. (38) 

Using 

m 2 = m § - e Am 2 tm - i Am|, g 
ml = m 2 -eAm 2 tm + |Am|, 
the inequality ( 1381) takes the form 

A; 4 m 4 + 2A; 2 m 2 + A; < 0, (40) 
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where 



/.'; 

k 2 

kn 



1 — 2(A + B) + (A — B) 
[A + B — (A — B) 2 ] eAr, 



atm 



(A - B) (1 - A-B) Ami 



2 



{A - B) eAml tm + -(A + B) Ami 



(41) 
(42) 

(43) 



Since k > 0, the inequality (|4"0"1) does not tolerate m 3 = 0; hence, there is a non-trivial 
lower bound on the neutrino masses. We want to find the numerical value of that bound. 
Using the values of A and B in equations (13"?]) . one finds that 



ki 

kn 



-16cV (1 



cos x) , 

2 (-2 + 13cV - 4cV cosx) eAm 2 tm 
+3 (c 2 - s 2 ) (-2 + 5cV + 4cV cos x) A 



in 



0' 



n 2 



3 (c 2 - s 2 ) eAm 2 tm + - (5 - lOcV - 8c 2 s 2 cosx) Am 2 Q 



(44) 

(45) 
(46) 



The case of real x, y, and w corresponds to cosx = +1- In (and only in) that case, the 
left-hand side of the inequality (1401) becomes linear in ra§; besides, in that case k 2 vanishes 
when c 2 s 2 = 2/9, thereby generating singularities at the points s 2 = 1/3 and s 2 = 2/3, as 
we saw in the previous section. 

For cosx 7^ +1, k A is negative. Since k is always positive, the inequality (j4"0|) then 
yields 



Vk 2 2 + 


ki 


k + k 2 







L. 



(47) 



The task now consists in finding the minimum value of L as a function of cos x (and of 
e = ±1); that minimum value provides the lower bound on m\. It is easy to convince 
oneself that L always has its minimum when cosx = — 1, for all experimentally allowed 

and AttIq. Computing L as a function of s 2 for fixed cosx = — 1, 
7.9 x 10~ 5 eV, we conclude the following: 



values of s 2 , Am 2 tm , 



Am atm 



2.5 x 1(T 3 eV, and Am 2 



• When the neutrino mass spectrum is normal, i.e. when e = +1, the minimum value 
of the lowest neutrino mass, mi, hardly varies with s 2 . One has mi > 1.679 x 10^ 2 eV 
for s 2 = 0.24 and m x > 1.665 x 10" 2 eV for s 2 = 0.40. 

• When the neutrino mass spectrum is inverted, i.e. when e = —1, the minimum 
value of the lowest neutrino mass, 777.3, varies strongly as a function of s 2 . One has 
m 3 > 2.9 x 10- 2 eV for s 2 = 0.24, m 3 > 9 x 10" 3 eV for s 2 = 0.40. 

Thus, the mass matrix of equation (fl~3|) with complex x, y, and w is not very predictive: 
it only allows one to derive a rather mild lower bound on the neutrino masses. There is 
also no prediction for the effective mass mpp, except for the rather trivial bounds 



1 2 2 

m 1 c 



m\s 2 \ < mpp < \m\c A + m A 2 s 



2„2 



,2„2| 



(4f 
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6 Conclusions 



In this paper we have constructed an extension of the Standard Model with three Higgs 
doublets 4>a an d four scalar gauge triplets A a and A 4 . The scalar triplets generate a 
type-II seesaw mechanism, thus explaining the smallness of the neutrino masses. We have 
employed a large horizontal symmetry group G, generated by the permutation group S3 
of the indices a and by six cyclic groups After spontaneous symmetry breaking, the 
charged-lepton mass matrix is diagonal; the different VEVs of the <p a allow for different 
charged- lepton masses m a (a = e, /x, r). In order to obtain a realistic neutrino mass 
matrix Ai u , we additionally allow for soft breaking of G, through terms of dimension two 
in the scalar potential. A crucial feature of our model is the equality among the diagonal 
entries of Ai u — this is one of the reasons for the predictiveness of the model. 

There are two relevant options: breaking G softly in the mass matrix of the scalar 
triplets either fully or keeping a /x <-> r symmetry intact; and having either hard or 
spontaneous CP breaking. Our model has the interesting property that spontaneous CP 
violation has no effect on M. U) i.e. it does not generate any physical phases in lepton mix- 
ing. The most predictive scenario combines the preservation of /x-r interchange symmetry 
with spontaneous CP violation, in which case we arrive at a viable neutrino mass matrix 
which has only three (real) parameters. This neutrino mass matrix leads to the usual pre- 
dictions of /x-r symmetric neutrino mass matrices, namely maximal atmospheric mixing 
and 9\z = — hence no CP violation in neutrino oscillations. Besides, the CP property 
mentioned before also prevents Majorana phases in our case. 

The solar mixing angle 9 is undetermined. Our three-parameter neutrino mass matrix 
predicts the neutrino masses rrij as functions of the two mass-squared differences and of 
6. For sin 2 9 < 1/3, which seems to be preferred by the data, we have a normal spectrum, 
while for sin 2 # > 1/3 the neutrino mass spectrum is inverted. When sin 2 9 — > 1/3 all 
the rrij diverge — see figure [TJ As for the effective mass mpp of neutrinoless 2f3 decay, our 
three-parameter mass matrix predicts mpp ~ m 3 /3. 

Acknowledgements: We thank Ernest Ma for suggestions and discussions that con- 
tributed significantly to this work. W.G. is grateful to H. Urbantke for discussions on finite 
groups. The work of L.L. was supported by the Portuguese Fundagao para a Ciencia e a 
Tecnologia through the projects POCTI/FNU/44409/2002, POCTI/FP/63415/2005, and 
U777-Plurianual. 

A The group structure of our model 

In this appendix we attempt a mathematical description of the full symmetry group of 
our model and of its irreducible representations (irreps). Clearly, S3 commutes neither 
with the z« of equation (j4]) nor with the Za of equation (jSJ), thus the full symmetry 
group is rather complicated. 
Let us define 

ni = diag(-l, 1, 1) , n 2 = diag(l, -1, 1) , n 3 = diag (1, 1, -1) , (Al) 
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(A2) 

h = ( 1 | , t 2 = ( 1 | , * 3 = ( 1 | . (A3) 
Then, 

N = {1, m, n 2 , n 3 , mn 2 , n 2 n 3 , n 3 ni, -1} (A4) 
forms an Abelian group isomorphic to Z 2 x Z 2 x Z 2 . Also, 

S 3 = {1, c+, c_, £i, t 2 , t 3 } (A5) 

forms a three-dimensional (reducible) representation of S3. 

Let us call G? the symmetry group utilized in this paper. G may be defined to be the 
group of the 6x6 matrices 



ms 
ns 



m,neN, seS 3 . (A6) 



This defining reducible representation of G may be called the 6. Clearly, G has 8x8x6 = 
384 elements^ Calling (m, n, s) the abstract element of G which is represented in the 6 
by the matrix of equation (IA6[) . the group multiplication law is 



(mi, m, si) (m 2 , n 2 , s 2 ) = (miSimaSj x , niSira 2 s 1 \ SaS 2 ) . (A7) 
From this group multiplication law it follows that G has eight one- dimensional irreps: 

l> 9 ' r ) : (to, n, s) -> (det m) p (det nf (det s) r , with p, q, r £ {0, 1 } . ( A8) 



It is obvious from equation (IA6I) that the matrices ms give a three-dimensional irrep of G, 
and similarly with the matrices ns. The matrices mns give one further three-dimensional 
irrep of G, since 

miniSim 2 n 2 s 2 = m 3 n 3 s 3 , with m 3 = miSim 2 s]~ 1 , n 3 = niSin 2 sf 1 , and s 3 = SiS 2 

(A9) 

complies with the multiplication law (1A7I) . Thus, G has 24 three-dimensional irreps: 



3i 



'' ' • (to, n, s) —y (det m) p (det n) q (det s) r ms, (A10) 
(m, n, s) — > (det m) p (det n) 9 (det s) r ns, (All) 
(m, n, s) — > (det m) p (det n) 9 (det s) r wins, (A12) 



o(p,q,r) 

H 2 



with p, q, r G {0, 1}. 



7 The 8 x 6 = 48 matrices ms, where m £ N and s S S3, form the Coxeter group B3. (We thank 
E. Ma for drawing our attention to Coxeter groups.) We may write G = N x B3 = (N x N) x 53, the 
symbol x denoting a semi-direct product. 
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The three-dimensional representations of G that we employ in our model are 

ms for (0 e , M , <j> T ) , 

runs for (e R , fi R , t r ) , 

ns for (D Le , D L/I , D Lt ) , 

(det n)ns for (A e , A M , A r ) . 

Our group G also has two and six-dimensional irreps, which are not used in our model. 
Next we include, for completeness, their construction. 

The group S3 has a two-dimensional irrep D 2 , generated by 

where u = (— 1 + / 2. Note that (det s) D 2 (s) is isomorphic to D 2 (s): 

(det 8)D 2 (s)= (J (A15) 

Therefore, G has four two-dimensional irreps: 

2 (p ' 9) : (m, n, s) -> (det m) p (det nf L> 2 (s) , with p, g G {0, 1} . (A16) 

The remaining irreps of G are four six-dimensional ones, 

6 M : (m,n,s) -> (det mf (det n) 9 D 6 (m, n, s) , with p,qE {0,1}- (A17) 

The irrep L>6 ( m ? ^, s ) is found in the decomposition of the product of the irreps ms and 
ns. Suppose there is a space <D 3 spanned by 61,2,3 transforming like ms, and another 
space <D 3 spanned by e[ 2 3 transforming like ns. Then, the space spanned by the <E> e' k 
(k = 1, 2, 3) transforms like mns, while the Cj <g> e' k with j 7^ k span a space which 
transforms like D 6 (to, n, s). It can be shown that this representation _D 6 (m, n, s) of G is 
irreducible, and also that it is equivalent to (det s) D 6 (to, n, s). 

The group G has the interesting property that it has no faithful irreps. It is obvious 
that the irreps with dimensions three and lower are not faithful. The six-dimensional 
irreps are not faithful either, as we now explain. Defining the elements a and b of G by 
a = (—1,1,1) and b = (1, —1, 1), then a, b, and ab generate the subgroups 1i 2 a \ 1> 2 \ 
and 1i 2 ab ^ of G, respectively. The isomorphisms 

6(0,0) ^ q 1 z (ab)^ g(i )0 ) ^ q 1 z (a) ) g( 0) l) ^ q j jifi) ^ 3(1,1) ^ Q j ^(a) y ^ 

(A18) 

are easy to demonstrate. Thus, none of the six-dimensional irreps represents G faithfully. 



B Spontaneous breaking of the fi— r symmetry 

Let us consider a simplified model with only two VEVs, t> M and v T . We assume the 
following symmetries: 

zi ■ -Vp, v T -> v T ; 

z 2 : Vp, v T -> -v T \ (Bl) 

Z 3 ■ V T . 
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The symmetries z\$ are assumed to be softly broken by terms of dimension two, while z% 
is assumed to be exactly conserved. For the sake of clarity we also assume all coefficients 
to be real. Then, 

V = a (k| 2 + H 2 ) + b (v;v T + + A (k| 2 + H 2 ) 2 + A' k| 2 H 2 . (B2) 

Only the b term breaks softly. 

Without loss of generality we take to be real and positive, writing 

— v cos0, v T = v sin <fi e ia , (B3) 

with v > and in the first quadrant. Then we obtain 

AV 4 

V = av 2 + bv 2 sin 20 cos a + Az/ 4 + — — sin 2 20. (B4) 



We require that 

dV 







5(20) 



/ AV 2 \ 

(V 2 cos 20) ( & cos cH — — sin 20 j . (B5) 



The solution cos 20 = corresponds to k| = \v T \ and is undesirable. But there is another 
solution, 

26 cos a 
^V 2 

which we adopt. Since the minimization of Vo in equation ( 1B4I) with respect to a leads 
to b cos a = — \b\ being negative, we must assume A' to be positive. If 



sm2( ^ = v^~. ( B6 ) 



|6| « (B7) 

which corresponds to the soft-breaking term being very small, then sin 20 <C 1 and \v^\ <C 
kl can be realized. 
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